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The length spectra ol flat three-dimensional dielectric resonators of circular shape were deter- 
mined from a microwave experiment. They were compared to a semiclassical trace formula obtained 
within a two-dimensional model based on the effective index of refraction approximation and a good 
agreement was found. It was necessary to take into account the dispersion of the effective index 
of refraction for the two-dimensional approximation. Furthermore, small deviations between the 
experimental length spectrum and the trace formula prediction were attributed to the systematic 
error of the effective index of refraction approximation. In summary, the methods developed in this 
article enable the application of the trace formula for two-dimensional dielectric resonators also to 
realistic, flat three-dimensional dielectric microcavities and -lasers, allowing for the interpretation 
of their spectra in terms of classical periodic orbits. 
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I. INTRODUCTION 

Open dielectric resonators have received great atten- 
tion due to numerous applications [1], e.g., as micro- 
lasers |2{5] or as sensors [BHH] , and as paradigms of open 
wave-chaotic systems [S^. The size of dielectric micro- 
cavities typically ranges from a few to several hundreds 
of wavelengths. Wave-dynamical systems that are large 
compared to the typical wavelength have been treated 
successfully with semiclassical methods. These provide 
approximate solutions in terms of properties of the cor- 
responding classical system. In the case of dielectric cavi- 
ties, the corresponding classical system is an open dielec- 
tric billiard. Inside the billiard rays travel freely while, 
when impinging the boundary, they are partially reflected 
and refracted to the outside according to Snell's law and 
the Fresnel formulas. The field distributions of resonance 
states of dielectric cavities can be localized on the pe- 
riodic orbits (POs) of the corresponding billiard [TOHTS] 
and the far-field characteristics of microlasers can be pre- 
dicted from its ray dynamics [I11I16) . Semiclassical cor- 
rections to the ray picture due to the Goos-Hanchen shift 
[TT] , Fresnel filtering [T8] , and curved boundaries [T9| are 
under investigation for a more precise understanding of 
the connections between ray and wave dynamics pUH^ . 
One of the most important tools in semiclassical physics 
are trace formulas, which relate the density of states of 
a quantum or wave-dynamical system to the POs of the 
corresponding classical system ^23r.25j . Recently, a trace 
formula for two-dimensional (2D) dielectric resonators 
was developed [26l [27] . The trace formula was success- 
fully tested for resonator shapes with regular classical 
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dynamics in experiments with 2D dielectric microwave 
resonators [55] and with polymer microlasers of various 
shapes |21]. However, typical microlasers like those 
used in Refs. [MS] are three-dimensional (3D) systems. 
While trace formulas for closed 3D electromagnetic res- 
onators have been derived [301 131] and tested [35] , hith- 
erto there is practically no investigation of the trace for- 
mula for 3D dielectric resonators. The main reason is the 
difficulty of the numerical solution of the full 3D Maxwell 
equations for real dielectric cavities. The case of flat 
microlasers is special since their in-plane extensions are 
large compared to the typical wavelength, whereas their 
height is smaller than or of the order of the wavelength. 
Even in this case complete numerical solutions are rarely 
performed. In practice, fiat dielectric cavities are treated 
as 2D systems by introducing a so-called effective index 
of refraction [SI (331 (see below) . This approximation has 
been used in Refs. [5, 29J and a good overall agreement 
between the experiments and the theory was found. How- 
ever, it is known [34] that this 2D approximation (called 
the rieff model in the following) introduces certain uncon- 
trolled errors. Even the separation between transverse 
electric and transverse magnetic polarizations intrinsic 
in this approach is not, strictly speaking, valid for 3D 
cavities |3S]. To the best of the authors' knowledge, no a 
priori estimates of such errors are known even when the 
cavity height is much smaller than the wavelength. The 
purpose of the present work is the careful comparison of 
the experimental length spectra and the trace formula 
computed within the 2D ricft approximation. Further- 
more, the effect of the dispersion of the effective index of 
refraction on the trace formula is investigated as well as 
the need for higher-order corrections of the trace formula 
due to, e.g., curvature effects. The experiments were per- 
formed with two dielectric microwave resonators of circu- 
lar shape and different thickness like in Ref. [M]. These 
are known to be ideal testbeds for the investigation of 
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wave-dynamical chaos [5S1[37] and have been used, e.g., 
in Refs. [28l [Ml [Ml [39]. The results of these microwave 
experiments can be directly applied to dielectric micro- 
cavities in the optical frequency regime if the ratio of the 
typical wavelength and the resonator extensions are sim- 
ilar. The paper is organized as follows. The experimen- 
tal setup and the measured frequency spectrum are dis- 
cussed in Sec. ini Section [lllj summarizes the UcB model 
for flat 3D resonators, the semiclassical trace formula for 
2D resonators and how these are combined here. The 
experimental length spectra are compared to this model 
in Sec. IIVI and Sec. \V\ concludes with a discussion of the 
results. 



frequency /. Plotting |S'2iP versus / yields the frequency 
spectrum. 

The measured frequency spectrum of disk B is shown in 
Fig. [2] It consists of several series of roughly equidistant 
resonances. The associated modes can be labeled with 
their polarization and the quantum numbers of the circle 
resonator, which are indicated in Fig. [2] Here, {m,nr) 
are the azimuthal and radial quantum number, respec- 
tively; TM denotes transverse magnetic polarization with 
the z component of the magnetic field, B^, equal to zero; 
and TE denotes transverse electric polarization with the 
z component of the electric field, E^, equal to zero. Each 
series of resonances consists of modes with the same po- 



ll. EXPERIMENTAL SETUP FOR THE 
MEASUREMENT OF FREQUENCY SPECTRA 

Two flat circular disks made of Teflon were used as 
microwave resonators. The first one, disk A, has a radius 
of i? = 274.8 mm and a thickness of 6 = 16.7 mm so 
R/b = 16.5. Its index of refraction is n = 1.434. A typ- 
ical frequency of / = 10 GHz corresponds to kR « 58, 
where k = 27r//c is the wave number and c the speed 
of light. The second one, disk B, has R — 225.1 mm, 
b = 10.1 mm {R/b = 22.3), and n = 1.430, with 10 GHz 
corresponding to kR w 47. The values of the indices of 
refraction n of both disks were measured independently 
(see Ref. [34]) and validated by numerical calculations 
[10] ■ They showed negligible dispe rsion in the consid- 
ered frequency range Figure 1(a) shows a photograph 
of the experimental setup. The circular Teflon disk is 
supported by three pillars arranged in a triangle. The 
prevalent modes observed experimentally are whispering 
gallery modes (WGMs) that are localized close to the 
boundary of the disk [34] . Therefore, the pillars perturb 
the resonator only negligibly because their position is far 
away from the boundary. Additionally, 4 cm of a spe- 
cial foam (Rohacell 31IG by Evonik Industries [H]) with 
an index of refraction of 71 « 1.02 is placed between the 
pillars and the disk as isolation [see Fig. 1(b)) . The to- 
tal height of the pillars is 260 mm so the resonator is 
not influenced by the optical table. Two vertical wire 
antennas are placed diametrically at the cylindrical side- 
walls of the disk (cf. Ref. [31]). They are connected to 
a vectorial network analyzer (PNA N5230A by Agilent 
Technologies) with coaxial rf cables. The network an- 
alyzer measures the complex scattering matrix element 
S'2i(/), where 



l^2l(/)|^ 



P2 



Pl,in 



(1) 



is the ratio between the powers P2,out coupled out via 
antenna 2 and Pi. in coupled in via antenna 1 for a given 



1 It was estimated that ^ X 10 GHz < 1.5 ■ 10 
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FIG. 1. (Color onli ne) |(a)| Photograph of the experimental 
setup with disk B. |(b)| Sketch of the setup (not to scale). 
Three metal pillars are used to support a Teflon disk with 
radius R and thickness b. A special foam (see text) is used 
to isolate the disk from the pillars. Two vertical wire an- 
tennas protruding from coaxial RF cables are used to couple 
microwave power into and out of the disk. A metal plate with 
variable distance D to the disk (not depicted in the photo) 
can be added to determine the polarization of the resonances. 



3 




12.0 12.2 12.4 12.6 12.8 13.0 24.5 24.7 24.9 25.1 25.3 25.5 

Frequency (GHz) 



FIG. 2. Measured frequency spectrum of disk B. The polarization and quantum numbers are indicated by TM or TE {m,nr). 
In the right part, only the relevant TM modes are indicated. The free spectral range (FSR) is the spacing between resonances 
with the same polarization and radial quantum number. 



larization and radial quantum number. The free spectral 
range (FSR) for each series is in the range of 145-170 
MHz. Only modes with rir < 5, that is, WGMs, are ob- 
served in the experiment. The quantum numbers were 
determined from the intensity distributions, which were 
measured with the perturbation body method (see Ref. 
[34j and references therein). To determine the polariza- 
tion of the modes a metal plate was placed parallel to the 
resonator at a variable distance D [see Fig. |l(b)| . Fig- 
ure [3] shows a part of the frequency spectrum with two 
resonances for different distances of the metal plate to 
the disk. The metal plate induces a shift of the reso- 
nance frequencies, where the magnitude of the frequency 
shift increases with decreasing distance D. Notably, the 
direction of the shift depends on the polarization of the 
corresponding mode: TE modes are shifted to higher fre- 
quencies and TM modes to lower ones, so the polarization 
of each mode can be determined uniquely. This behavior 
is attributed to the different boundary conditions for the 
field (TM modes) and the field (TE modes) at 
the metal plates. The former obeys Dirichlet, the latter 
Neumann boundary conditions. A detailed explanation 
is given in Appendix [A) 



0.08 



0.06- 



^ 0.04- 



0.02- 




12.20 



12.22 12.24 12.26 
Frequency (GHz) 



12.28 12.30 



FIG. 3. Frequency spectrum without metal plate (solid line) 
and with metal plate [see Fig. 1(b) at distances D = 14 mm 
(dashed), 10 mm (dotted), and 8 mm (dot-dashed) for disk 
B. With decreasing distance D the TE modes are shifted to 
higher frequencies, while the TM modes are shifted to lower 
frequencies. 



III. THE EFFECTIVE INDEX OF 
REFRACTION AND THE TRACE FORMULA 

The open dielectric resonators are described by the vec- 
torial Helmholtz equation 



= 



(2) 



with outgoing-wave boundary conditions, where E and B 
are the electric and the magnetic field, respectively, and 
n{r) is the index of refraction at the position f. Though 



all components of the electric and magnetic fields obey 
the same Helmholtz equation, they are not independent 
but rather coupled in the bulk and at the boundaries as 
required by the Maxwell equations. The eigenvalues kj of 
Eq. ([2| are complex and the real part of kj corresponds 
to the resonance frequency fj = c Re {kj ) / {2tt) of the 
resonance j, while the imaginary part corresponds to the 
resonance width Tj = —clm{kj) /tt (full width at half 
maximum) . 

For the infinite slab geometry (see Fig. |4|, the vectorial 
Helmholtz equation can be simplified by separating the 
wave vector k into a vertical z component, kze^, and a 
component parallel to the x-y plane, fcn . Thus, k^ = {k'^ + 
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FIG. 4. Ray traveling through an infinite dielectric slab with 
thickness 6 and index of refraction n. The wave vector k 
is decomposed into its components perpendicular (kzCz) and 
parallel (fc||) to the plane of the disk, where fc|| = riesk. The 
angle of incidence on the top and bottom surface is 9. 

/C| and the angle of incidence on the top and bottom 
surface of the resonator \s = arctan (fcy/Zc^). For a 
resonant wave inside the slab the wave vector component 
kz must obey the resonance condition 



exp(2ifc26)r2(6') = 1^ 



(3) 



where r is the Fresnel reflection coefficient. The solutions 
of Eq. (|3| yield the quantized values of kz ■ The effective 
index of refraction riofi is deflned as 



k\\ 



(4) 



and corresponds to the phase velocity with respect to the 
x-y plane. In the experiments only modes trapped due 
to total internal reflection (TIR) are observed. In this 
case, the reflection coefflcient r can be written as 



r = exp (— 2i(5o) 



with 



tan^o = vi. 



(5) 



(6) 



'off 



where v = for TM modes and v = 1 for TE 
modes. With these definitions, the quantization condi- 
tion Eq. (|3| for k^ can be reformulated as an implicit 
equation for the determination of Uoff , 
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"off 



2 arctan v 




(7) 



with C = 0, 1, 2, . . . being the order of excitation in the 
z direction [S]. The arctan term in Eq. ([?]) corresponds 
to the Fresnel phase due to the reflections and the Qti 
term to the geometrical phase. In the framework of the 
neft model the flat resonator is treated as a dielectric slab 
waveguide and the vectorial Helmholtz equation [Eq. ^] 
is accordingly reduced to the 2D scalar Helmholtz equa- 
tion [3 [331 IMj by replacing n by the effective index of 
refraction Uoff, 
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FIG. 5. Effective index of refraction for the TM modes with 
^ = of disk A (solid line) and disk B (dashed line). 



where the wave function ^in,out inside, respectively, out- 
side of the resonator corresponds to Ez in the case of TM 
modes and to Bz in the case of TE modes. The bound- 
ary conditions at the boundary of the resonator in the 
x-y plane (i.e., the cylindrical sidewalls), dS, are 



4-; 



in|9S 



out \ dS 



and 



dfi 



as 



df 



(9) 



dS 



where n is the unit normal vector for dS, ^ — I for 
TM modes, and /i = i/ri^ff for TE modes. Equation 
([s]) can be solved analytically for a circular dielectric res- 
onator [42 . However, it should be stressed that Eq. ([s]) 
is not exact for flat 3D cavities. It defines the 2D nog 
approximation whose accuracy is unknown analytically 
but which has been determined experimentally in Ref. 
j34j . Our purpose is to investigate the precision of this 
approximation for the length spectra of simple 3D dielec- 
tric cavities. The effective index of refraction for the TM 
modes with the lowest z excitation C = of disk A and 
B is shown in Fig. [5] Obviously, ricff depends strongly 
on the frequency, and this dispersion plays a crucial role 
in the present work. It should be noted that also TE 
modes and modes with higher z excitation exist in the 
considered frequency range, however, in the following we 
focus on TM modes. 

The density of states (DOS) in a dielectric resonator 
is given by [26] 



p{k) 



TT ^ [A: - Re (kj)]^ + [Im (kj)]"^ 



(10) 



where the summation runs over all resonances j. The 
DOS can be separated into a smooth, average part p 
and a fiuctuating part, pfluc- The smooth part is well 
described by the Weyl formula given in Ref. [221 and de- 
pends only on the area, the circumference and the index 
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of refraction of the resonator. The fluctuating part, on 
the other hand, is related to the POs of the corresponding 
classical dielectric billiard. For a 2D dielectric resonator 
with regular classical dynamics, the semiclassical appro- 
ximation for pfluc is [26j 



Pfluc(fc) 




Spo|i?po|V^e'(""-+^-)+c.c. (11) 



where ipo is the length of the PO, Rpo is the product of 
the Fresnel reflection coefficients for the reflections of the 
rays at the dielectric interfaces, i^po denotes the phase 
changes accumulated at the reflections [i.e., arg(i?po)] 
and at the conjugate points of the corresponding PO, and 
the amplitude Bpo is proportional to Apo/ y^lpo, where 
Apo is the area of the billiard covered by the family of 
the PO. It should be noted that this semiclassical for- 
mula fails to accurately describe contributions of POs 
with angle of incidence close to the critical angle for TIR, 
ttcrit = arcsin(l/n), as concerns the amplitude. Con- 
sequently, higher-order corrections to the trace formula 
need to be developed for these cases [55]. We restrict 
the discussion to the experimentally observed TM modes 
and compare the results with the trace formula obtained 
in Ref. W . To select the TM modes with C = the 
polarization and the quantum numbers of all resonan- 
ces had to be determined experimentally as described in 
Sec. ini The trace formula for 2D resonators is applied 
to the flat 3D resonators considered here by inserting the 
frequency-dependent effective index of refraction ncs{k) 
instead of n into Eq. (111. To test the accuracy of the 



resulting trace formula we computed the Fourier trans- 
form (FT) of pfluc- In Ref. [5] it was shown that it is 
essential to fully take into account the dispersion of rias 
in the FT for a meaningful comparison of the resulting 
length spectrum with the geometric lengths of the POs. 
Therefore, we define 



P(^) = / dfcpfluc(A:)exp{-zfcncff(fc)^} 

fe„i„ (12) 
= j:exp{^tkjn,s{kj)e} -FT{p{k)} , 

3 

where the quantity ^ is a geometrical length and \p{() \ is, 
thus, called the length spectrum. We will compare it to 
the FT as defined in Eq. ( |l2| ) using nes{k) instead of n 
of the trace formula, /5sci(^)- The resonance parameters 
kj are obtained by fitting Lorentzians to the measured 



frequency spectra, and /cmin.max = 27r/i„in,max/c corre- 
spond to the frequency range considered. Since in a circle 
resonator the resonance modes with m > are doubly 
degenerate, the measured resonances are counted twice 
each. Note that even though only the most long-lived re- 
sonances (i.e., the WGMs) are observed experimentally, 
and these comprise only a fraction of all resonances, a 
comparison of the experimental length spectrum with the 
trace formula is meaningful |28| . 



IV. COMPARISON OF EXPERIMENTAL 
LENGTH SPECTRA AND TRACE FORMULA 
PREDICTIONS 



Figure [6] shows the experimental length spectrum eval- 
uated using Eq. ( 12 1 and the FT of the semiclassical trace 
formula, Eq. (11), for disk A. A total of 572 measured 
TM modes with radial quantum numbers rir = 1-5 from 
/min = 6.8 GHz to /max = 20.0 GHz was used. The POs 



in the circle billiard are denoted by their period q and 
their rotation number 77 and have the shape of polygons 
and stars (see insets in Fig.[6|. Their lengths are indi- 
cated by the solid arrows. The POs with 2 < g < 15 and 
77 = 1 were used to compute the trace formula. Only POs 
with q < 8 are indicated in Fig. [6] for the sake of clarity. 
The POs with q > 9 only add small contributions to the 
right shoulder of the peak corresponding to the (7, 1) and 
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FIG. 6. Length spectrum for the TM modes of disk A. The 
solid line in the top part is the experimental length spectrum 
and the dashed line the FT of the semiclassical trace formula. 
The curves in the bottom part are the contributions of the 
individual POs to the semiclassical trace formula. The solid 
arrows indicate the lengths of the POs denoted by {q, rj) and 
the dash ed arrows indicate the peak positions computed from 
Eq. | |17[ ). The circumference 2ttR of the circle is also indi- 
cated. 



6 



90 
80 
70 



^ 60 



■g 50 

e 

40 
30 
20, 





1 1 


\ ^ 

; V 

V (6,1) 
(5_i) 

; ! 


. apo(5, 1) ; 




: 1 ! T— 1 


(4, 1) 


U| i 4 


(3,1) 







5 10 15 

Frequency (GHz) 



20 



FIG. 7. Critical angle Ocrit = arcsin (l/ricft) for TIR with re- 
spect to the frequency for the TM modes of disk A. The angles 
of incidence of the {q, rj) orbits are indicated by the horizontal 
lines, and the critical frequencies fcrit at which they become 
confined by TIR by the vertical lines. The (3, 1) and (4, 1) or- 
bits are not confined for the frequency range considered here. 



1) orbits. The amplitudes are 



= 7^"^^^° = i 2; othe;;ise ' (13) 



1, q = 2Tj 



and the phases are 



t^po = ^ - g 2 + 9 ['■("po) 



(14) 



with the angles of incidence of the POs (with respect to 
the surface normal) being 



7r/2 — rjTr/q . 



(15) 



The overall agreement between the experimental length 
spectrum and the semiclassical trace formula is good and 
the major peaks in the length spectrum are close to the 
lengths of the (5, l)-(8, 1) orbits. However, no clear peaks 
are observed at the lengths of the (3, 1) and the (4, 1) or- 
bit in the experimental length spectrum. Note that in the 
experimental length spectrum only orbits that are con- 
fined by TIR are observed (cf. Ref. [55] )• This is not the 
case for the (3, 1) and the (4, 1) orbits in the frequency 
range of interest, where their angle of incidence apo is 
smaller than the critical angle acrit = arcsin (l/ricff) as 
depicted in Fig. [7j The length spectrum of disk B is 
shown in Fig. [8j Altogether 546 resonances with = 1- 
5 from /mill = 9.8 GHz to /max = 26.9 GHz were used. 
The semiclassical trace formula was again computed for 
2 < q < 15 and 77 = 1. The agreement of the experi- 
mental length spectrum and the FT of the trace formula 
is good and comparable to that obtained for disk A. As 
in the case of disk A, the experimental length spectrum 
exhibits no peaks for the (3, 1) orbit, whose length is not 
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FIG. 8. Length spectrum for the TM modes of disk B. The 
solid line in the top part is the experimental length spectrum 
and the dashed line the FT of the semiclassical trace formula. 
The curves in the bottom part are the contributions of the 
individual POs to the semiclassical trace formula. The solid 
arrows indicate the lengths of the POs denoted by {q, if) and 
the dashed arrows indicate the peak positions computed from 
Eq. | |17[ ). The circumference 211 R of the circle is also indi- 
cated. 



within the range depicted in Fig. [8) and for the (4, 1) 
orbit since they are not confined by TIR in the consid- 
ered frequency range. A closer inspection of Figs. [6] and 
|8] shows two unexpected effects. First, the peak posi- 
tions of the FT of the semiclassical trace formula deviate 
slightly from the lengths of the POs. This can be seen 
best in the bottom parts of Figs. [6]and|8j where the con- 
tributions of the individual POs to the trace formula are 
depicted. Second, there is a small but systematic dif- 
ference between the peak positions of the experimental 
length spectrum and those of the FT of the trace formula. 
We will demonstrate that the first effect is related to the 
dispersion of ricft and the second effect to the systematic 
error of the ncff model. 

The difference between the peak positions of the trace 
formula and the lengths of the POs can be understood 
by considering the exponential term in the FT of the 
semiclassical trace formula Eq. (Ill, which for a single 
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TABLE I. Comparison of the PO lengths ipo, the estimates for 
the peak positions £pcak from Eq. I\17\ , and the peak positions 
^pcak of the single orbit contributions to the semiclassical trace 
formula (see bottom parts of Figs. [6| and [s]) . The (9, 1) and 
the (10, 1) orbit are not indicated in Figs. |6] and [s] since they 
give only small contributions to the rightmost peak in each 
figure. The top part is for disk A, the bottom one for disk B. 
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PO is 



exp {i[ncs{k)k{£po - £) + ^po]} 



(16) 



with ippo given by Eq. ( |14[ ) . The crucial point is that the 
phase (fipo is frequency dependent because it contains the 
phase of the Fresnel coefhcients, which, in turn, depends 
on noff(fc). The modulus of the FT will be largest for that 
length i for which the exponent in Eq. ( 16 ) is stationary, 
i.e., its derivative with respect to k vanishes. This leads 
to the following estimate f peak for the peak position, 



^pcak — ^po '^Q 



dn, 



dk 



dk 



ko 



where 



(5(a,7ioff) — arctan ■ 



' n^g sin^ a — 1 
rics cos a 



(17) 



(18) 



is related to the Fresnel coefficients via Eq. ([5|. The wave 
number fcp at which Eq. (17) is evaluated is the center of 



the relevant wave number / frequency interval, which is 




^ciit ^inin 
^crit ^ ^min 



(19) 



with fccrit — 27r/crit/c. The frequency fait is defined by 
sinQ!po = l/'T-eff (/crit), i-G-, it corrcsponds to the mini- 
mum frequency at which the considered PO is confined by 



TIR (cf. Fig. [7]). Below fccrit the Fresnel phase vanishes. 
The estimated peak positions ^peak are indicated by the 
dashed arrows in Figs. [6] and |8] and agree well with the 
peak positions i'^^^]^ of the individual PO contributions in 
the bottom parts of the figures. A list of the lengths of the 
POs, the peak positions of the single PO contributions, 
and the estimates ^poak according to Eq. (17 1 for disks 
A and B is provided in Table ID In general, the estimate 
^pcak deviates only by 1-2 mm from the actual peak po- 



sition il^^^^y^ (about l%c 



oil 



scl 
peak 



). The (3,1) and the (4,1) 



orbits are not confined by TIR. Therefore, for these orbits 
the Fresnel phase vanishes and accordingly ^peak = ^po- 
Furthermore, their contributions to the length spectrum 
are symmetric with respect to ^ = ipo, while those of 
the other POs are asymmetric with an oscillating tail to 
the left (see bottom parts of Figs. [6]and[8]). These tails 
are attributed to the frequency dependence of the Fresnel 
phase. They can lead to interference effects, as can be 
seen for example in Fig. [8j There, e.g., the peak positions 
of the semiclassical trace formula (dashed line in the top 
part) for the (5, 1) and the (6, 1) orbit deviate from the 
peak positions ^pg^k of the corresponding single orbit con- 
tributions (dashed and dotted lines in the bottom part) 
due to interferences between the contributions of a PO 
and the side lobes of those of the other POs. In order to 
identify such interferences, it is generally instructive to 
compare the FT of the semiclassical trace formula with 
those of its single orbit contributions. It should be noted 
that the effect discussed in this paragraph also occurs for 
any 2D resonator made of a dispersive material. 

It was shown in the previous paragraph that the dis- 
persion of Tioff plays an important role. Furthermore, the 
semiclassical trace formula is known to be imprecise for 
POs with angles of incidence close to the critical angle. 
This is especially crucial here since several POs are close 
to the critical angle in at least a part of the considered 
frequency regime (see Fig. [7]). These deficiencies of the 
semiclassical trace formula indicate the necessity to im- 
plement modifications of it. To pursue this presumption 
we will compare the experimental length spectrum with 
the FT of the exact trace formula for the 2D dielectric cir- 
cle resonator using a frequency-dependent index of refrac- 
tion n(fc) in order to investigate the deviations between 
it and the FT of the semiclassical trace formula. The 
trace formula is called exact since it is derived directly 
from the quantization condition for the dielectric circle 
resonator and without semiclassical approximations. It 
is given by 



Pflu 



— OO 

(20) 



po=(g,i)) 
with the definitions 



E„i{x) 



H 



(1) 



H 



(2; 



■{nx) 



(21) 
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and 



Rm{a 



Bmix) 



h: 



'(1) 



H 



(1) 



(nx) 



h: 



'(1) 



H 



(1) 



(x) 



h: 

t — 

H 



'(2) 

m 

W) 



(nx) 



'(1) 

m 



(22) 



(23) 



(24) 



(25) 
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f 1 21 

Here, x = kR, Hm {x) are the Hankel functions of the 
first and second kinds, respectively, and the prime de- 
notes the derivative with respect to the argument. Equa- 
tion (20) is essentiaUy Eq. (67) of Ref. 26 with an ad- 



ditional factor (1 -I- ^^) in the term P„j(a;). A detailed 
derivation is given in Appendix [B] In the semiclassical 
limit, the term Wj^^ in Eq. (20) turns into the product of 



the Fresnel reflection coefficients, the term e^'^"'™£'^ into 
the oscillating term e^^-'^p", and Pm contributes to the 
amplitude Bpo- For POs close to the critical angle, i.e., 
when the stationary point of the integrand in Eq. ( 20 ) 
is mo ~ nfci? sin ttcrit , the term Rf^ varies rapidly with 
m, whereas it is assumed to change slowly in the station- 
ary phase approximation used to derive the semiclassical 
trace formula 26^. Therefore, including curvature cor- 
rections in the Fresnel coefficients does not suffice for an 
accurate calculation of the contributions of these POs to 
Pfluc- Consequently, we compute the integral entering in 
Eq. ( |20| ) numerically. 

Figure |9] shows the comparison of the experimental 
length spectrum (solid line), the FT of the semiclassi- 
cal trace formula (dashed line), Eq. ( 11 ), and that of the 
exact trace formula (dotted line), Eq. (20), for disk A. 
The latter was computed for 2 < q < 10 and rj = 1. 
The other two curves are the same as in the top part of 
Fig. [6] The main difference between the semiclassical and 
the exact trace formula are the larger peak amplitudes of 
the exact trace formula, where this is mainly attributed 
to the additional (1 -I- ^^) factor. The differences be- 
tween the peak amplitudes of the experimental length 
spectrum and those of the exact trace formula are ac- 
tually expected since the measured resonances comprise 
only a part of the whole spectrum (cf. Ref. [IH]); though 
they are not very large. On the other hand, the peak po- 
sitions ^max* of the exact trace formula differ only slightly 
from those of the semiclassical one, and still deviate from 
those of the experimental length spectrum, £™pJ (see in- 
set of Fig. [9]). The difference A^^ax = ^^ri - ^max* is in 
the range of 3-4 mm, i.e., about 2%c of the periodic orbit 
length £po. Similar effects are observed in Fig. [TO) which 
shows the experimental length spectrum and the FT of 
the semiclassical trace formula and that of the exact trace 



FIG. 9. Length spectrum for the TM modes of disk A. The 
solid line is the experimental length spectrum, the dashed 
line the FT of the semiclassical, and the dotted line the FT 
of the exact trace formula. The dashed arrows indicate the 
lengths £poak . A magnification around the peak corresponding 
to the (6, 1) orbit is shown in the inset. Here, A^max denotes 
the difference between the positions of the maxima of the 
experimental length spectrum and those of the FT of the exact 
trace formula. 
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FIG. 10. Length spectrum for the TM modes of disk B. The 
solid line is the experimental length spectrum, the dashed line 
the FT of the semiclassical, and the dotted line the FT of the 
exact trace formula. The dashed arrows indicate the lengths 
^pcak. The solid and dashed lines are the same as in the top 
part of Fig. [8] 



formula (computed for 2 < g < 10 and rj — 1) for disk 
B. The peak amplitudes of the exact trace formula are, 
again, somewhat larger than those of the experimental 
length spectrum, and the peak positions -^max' of the ex- 
act trace formula differ from those of the experimental 
length spectrum by A^max ~ 2 mm for the (5,1) and the 
(6, 1) orbits. The relative error ^^max/^max ~ 1.5%o is, 
thus, slightly smaller than in the case of disk A. 
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FIG. 11. (a) Difference between the calculated and the mea- 
sured resonance frequencies with respect to the measured fre- 
quency for the TM modes with rir = 1 of disk A. The modes 
with Ur > 1 show a similar behavior (adapted from Ref. [34]). 
(b) Measured (dots) and calculated (solid line) free spectral 
ranges for TM modes with = 1. The calculated FSR is 
plotted as a curve instead of data points to guide the eye. 
The inset shows a magnification in the frequency range of 
10-18 GHz. 



Since the trace formula itself is exact, the only explana- 
tion for these deviations is the systematic error of the rics 
model. Therefore, we compare the measured (/cxpt) a-nd 
the calculated (/caic) resonance frequencies for disk A in 
Fig. [TT] The resonance frequencies were calculated by 
solving the 2D Helmholtz equation [Eq. ([8|] for the cir- 
cle as in Ref. [31|. The difference between the measured 
and the calculated frequencies in Fig. [TTJa) is as large 
as half a FSR (i.e., the distance between two resonances 
with the same radial quantum number) and slowly de- 
creases with increasing frequency. This is in accordance 
with the result that the FSR of the calculated resonances 
in Fig. [Tl^b) is slightly larger than that of the measured 
ones. Since the frequency spectrum consists of series of 
almost equidistant resonances, the effect of this system- 
atic error on the peak positions in the length spectrum 
can be estimated by considering a simple ID system with 
equidistant resonances whose distance equals FSR. The 
peak position in the corresponding length spectrum is 
^max 1/FSR, and a deviation of AFSR leads to an 



AFSR 
FSR ' 



(26) 



of the peak position. With AFSR = FSR, 



'CXpt " 



-FSR, 



■calc 



-0.4 MHz compared to FSR « 120 MHz we expect a de- 
viation of 3%c or A^niax ~ 5 mm in the peak positions, 
which agrees quite well with the magnitude of the devia- 
tions found in Fig. [9] For disk B, the comparison between 
the measured and the calculated FSR (not shown here) 
yields AFSR/FSR « -l.l%c, which also agrees weU with 
the deviations of A^max/^max ~ 1.5%o found in Fig. \T0\ 
Thus, we may conclude that the deviations between the 
peak positions of the experimental length spectrum and 
that of the trace formula indeed arise from the system- 
atic error of the rioff model. Unfortunately, we know of no 
general method to estimate the magnitude of this system- 
atic error beforehand. It should be noted that the exact 
magnitude of the systematic error contributing to the de- 
viations found in Figs. [9| and [T0| depends on the index of 
refraction used in the calculations. Still, it was shown 
in Ref. [31] and also checked here that deviations remain 
regardless of the value of n used, which is known with per 
mill precision for the disks A and B piP. Furthermore, 
Fig. [Tl^b) demonstrates that the index of refraction of a 
disk cannot be determined without systematic error from 
the measured FSR even if the dispersion of rics is fully 
taken into account. 



V. CONCLUSIONS 

The resonance spectra of two circular dielectric mi- 
crowave resonators were measured and the corresponding 
length spectra were investigated. In contrast to previous 
experiments with 2D resonators [28], flat 3D resonators 
were used. The length spectra were compared to a combi- 
nation of the semiclassical trace formula for 2D dielectric 
resonators proposed in Ref. [26 and a 2D approximation 
of the Helmholtz equation for flat 3D resonators using an 
effective index of refraction (in accordance with Ref. [5 ) . 
The experimental length spectra and the trace formula 
showed good agreement, and the different contributions 
of the POs to the length spectra could be successfully 
identified. The positions of the peaks in the experimen- 
tal length spectrum are, however, slightly shifted with 
respect to the geometrical lengths of the POs. We found 
that this shift is related to two different effects, which 
are, first, the frequency dependence of the effective index 
of refraction and, second, a systematic inaccuracy of the 
TT-cff approximation. In the examples considered here, the 
former effect is as large as 5%c of the PO length while the 
latter effect is as large as 2%c of ^po, and the two effects 
cancel each other in part. The results and methods pre- 
sented here provide a refinement of the techniques used 
in Refs. [3 [29] and allow for the detailed understand- 
ing of the spectra of realistic microcavities and -lasers in 
terms of the 2D trace formula. Furthermore, many of 
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the effects discussed here also apply to 2D systems made 
of a dispersive material. Some open problems remain, 
though. The comparison of the semiclassical trace for- 
mula with the exact one for the circle showed that the 
former needs to be improved for POs close to the critical 
angle. Furthermore, there are some deviations between 
the experimental length spectra and the trace formula 
predictions due to the systematic error of the effective in- 
dex of refraction approximation. Its effect on the length 
spectra proved to be rather small and, thus, allowed for 
the identification of the different PO contributions. How- 
ever, the computation of the resonance frequencies of flat 
3D resonators based on the combination of the 2D trace 
formula and the rics approximation would lead to the 
same systematic deviations from the measured ones as 
in Ref. |34j . Another problem with this systematic error 
is that its magnitude cannot be estimated a priori. The 
comparison of the results for disk A and disk B seems to 
indicate that it gets smaller for b/R — > 0, but there are 
not enough data to draw final conclusions yet, especially 
since the value of b/R is of similar magnitude for both 
disks. In fact, Ref. [31] rather indicates that the sys- 
tematic error of the rics model increases with decreasing 
b/R. This could be attributed to diffraction effects at 
the boundary of the disks that become more important 
when b gets smaller compared to the wavelength. On the 
other hand, the exact 2D case is recovered for b — ^ oo, 
i.e., an infinitely long cylinder. In conclusion the accu- 
racy of the rtoff model in the limit b/R — > remains an 
open problem. An analytical approach to the problem 
of flat dielectric cavities that is more accurate than the 
ries model would be of great interest. Another perspec- 
tive direction is to consider 3D dielectric cavities with 
the size of all sides having the same order of magnitude 
and to develop a trace formula for them similar to those 
for metallic 3D cavities l30H32l. 
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Appendix A: Influence of the metallic plate on the 
resonance frequencies 

The placement of a metal plate parallel to the dielec- 
tric disk influences the effective index of refraction, thus, 
leading to a shift of the resonance frequencies. In order 
to determine the change of n^g, we, flrst, calculate the 
reflection coefficient r{D) for a wave traveling inside the 
dielectric and being reflected at the dielectric-air inter- 
face with the metal plate at a distance D. The geometry 



used here is depicted in Fig. [12] The ansatz for the Ez 
field (TM polarization) and, respectively, the field 
(TE polarization) is 

A«'(x,?/)(e^'''-^ -f re-*'=-^)e-*"* : z < 



0<z<D, 



where * fulfiUs (A -I- kp^ = 0, w = 2nf is the angular 
frequency. A, bi 2 are constants, and r is the reflection 
coefficient. The different wave numbers are connected by 



(A2) 



For the case of TIR considered here, q is real, and the 
penetration depths I of the field intensity into region II is 
I = 1/(2(7). For the case of TM polarization the electric 
field is 



= B* cosh [q{z - D)]e- 



(A3) 



since it obeys Neumann boundary conditions at the metal 
plate, where B is a constant. The boundary conditions 
at the interface between region I and II are that n?Ez 
and ^^1^ are continuous, which yields 

kz — in^gtanh (qD) 



kz + m^gtanh (qD) 



(A4) 



For the case of TE polarization the magnetic field in re- 
gion II is 



=B«'sinh [q{z-D)]e' 



(A5) 



since it obeys Dirichlet boundary conditions at the metal 
plate. With the condition that Bz and are continu- 
ous at the dielectric interface, 

kz — iq coth {qD) 



kz + igcoth (qD) 



(A6) 



is obtained. Analogous to Eq. (|6| we define r{D) 
exp [—2i5{D)] and obtain 



tan [d{D)] = v\ 



'cff 



-h[D/i2l)] (A7) 



Vff 



D 




metal 



II 



air, n = 1 




FIG. 12. Geometry and coordinate system for the calculation 
of the reflection coefScient r{D). 
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with 

^ I tanhx : TM polarization (A8) 
1 coth X : TE polarization . 

For D ^ I, Eq. ([6| is recovered. This explains why the 
optical table does not disturb the resonator. By inserting 
r{D) into Eq. ^ we obtain 



kb 



1 



-.{do + d{D) + Cn}^gin,s,D) (A9) 



''eff 



as condition for ncs{D), where So refers to Eq. ([6|. In 
order to determine the effect of a change of D on the 
effective index of refraction, we compute 



ducff 
dD 



1 d5{D) / dg 



^n^-^^ dD 



dn, 



eff 



where 



tan(5o 
2l[l + h^ tan2 (,5o); 

The derivative of h is 



dSjP) 
dD 



(AlO) 



(All) 



h'{x) 



cosh ^ (x) 



TM 



sinh~^{a;) : TE 



(A12) 



which approaches h'{x) = Ave for a; — >■ oo with v = 
+1 for TM and i; = — 1 for TE polarization. Then, for 



ducff 



C 



dD =-r' 



-D/l 



(A13) 



with C > given as 
C ' ' 



tan(5o 



dg 

9"ott 



v/n2 - 1 + tan2 {So) 



(A14) 



Accordingly, for large distances of the metal plate, 

' +Ce-^/^ : TM 



-Ce-^/' : TE, 



(A15) 



i.e., for TM modes rieff increases for decreasing distance 
D of the metal plate from the resonator and for TE 
modes it gets smaller. This qualitative behavior is found 
for the whole range of D, even when Eq. (A15) is no 



longer valid. Since the resonance frequencies of the 
disk are to first-order approximation c>c l/riog, the TM 
modes are shifted to lower and the TE modes to higher 
frequencies as observed in Fig. |3] 



Appendix B: Exact trace formula for the circular 
dielectric resonator 



The DOS for the TM modes of the 2D circular dielec- 
tric resonator is 



m— — oo n^ — l ^ 



(Bl) 



where the eigenvalues km. are the roots of 02] 

Sraix) = x[nyra{nx)Ri\\x) ~ J„, (7ix)H'|^^ (a:)] (B2) 

with X — kR. The factor x ensures that s„j has no poles. 
Informally one can write Smix) — /(x)n„^(a; — Rkm,nr)j 
where f{x) is a certain smooth function without zeros 
and poles. This means that the DOS can be written as 



p{k) 



p +00 

^ y Im 



dsm 
dx 



smooth term. (B3) 



m— — oo 



The smooth term contributes to the Weyl expansion and 
requires a separate treatment ^26l [S^ . In the following 
we ignore all such terms. The derivative in Eq. (B3) 
contains two terms. 



ds 
dx 



ds„ 



dx 



dsjn dn 
dn dx 



The first term is 
dsm 



dx 



= -x{n^ -l)im{nx)}i^^\x) 



(B4) 



(B5) 



where the second derivatives of the Bessel and Hankel 
functions were resolved via the Bessel differential equa- 
tion. The second term is 



n2 ( jjU) ) (^) 



(B6) 



where j^{nx) was replaced by ~ ^^i) (^) since Sm{x) = 0. 
We approximate 



h: 



'(1) 



h: 



(iy(^)«- 



- 1 



(B7) 



and obtain 



dsr. 



dn 



X 

n 



{n^ -l)3m{nx)Rl^\x). (B8) 



It was checked numerically that this approximation is 
very precise, which is why we still call the result an exact 
trace formula. Combining both terms yields 



dx 



{n^ - 1) 



k dn 
n dk 



(B9) 
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with 



H 



'(1) 



Sm(x) = nj^{nx) r[)(^) 

The DOS is, thus, 



(BIO) 



p{k) 



2m 



We replace the Bessel functions in by (Hm^ + 

Hm ) /2 and extract a term [Em {x) + 1] with Em defined 
in Eq. (21 1 to obtain 



E„ 



1 



[1 — EmRm] , 



(B12) 



where Rm, Am, and Bm are defined in Eqs. (23), (24 1, 



and (|25|), respectively. With the help of the geometric 

(B13) 



series l/Sm is rewritten as 

1 1 . °° 

^3 "^^ ^" Pm ^ ^ ( -^m Rm ) 



with 



9=1 



1 + — I s„ 



Using the Wronskian 



Ai 



TTZ 



(B14) 



(B15) 



this simplifies to 

p — 



4i 



(B16) 



and we obtain 



dk J 



+ 00 



q=l 



(B17) 

The first term, 1/ Bm, corresponds to the smooth part of 
the DOS. Since we are only interested in the fluctuating 
part, we drop the term and apply the Poisson resumma- 
tion formula to the rest to obtain 



+00 



?7— — oo_ 



+00 



/ dme^-'""WmY,{EmRnrY+C.C. 

(B18) 

+00 oo 

with Pm defined in Eq. ( [22| ). Replacing ^ ^ with 

r/— — cxD q—1 

2 ^ and ignoring those (g, 77) combinations that are not 

po 

related to any POs and, thus, do not give significant con- 
tributions finally yields Eq. ([20|). Taking the semiclassi- 



cal limit as described in Ref. |26 results in Eq. ( 11 ) with 
an additional factor of (1 + ^^). This means that the 
dispersion of n leads to slightly higher amplitudes in the 
semiclassical limit. 
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